arXiv:1505.00804v2 [cond-mat.str-el] 21 Dec 2015 


Superconducting gap and vortex lattice of the heavy fermion compound CeCu2Si2 


Mostafa Enayat,^’0 Zhixiang Sun,^’0 Ana Maldonado,^’ ^’0 Hermann Suderow,^’^ 
Silvia Seiro,^ Christoph Geibel,^ Steffen Wirth,^ Frank Steglich,^ and Peter Wahl^A,0 

^Max-Planck-Institut fiir Festkorperforschung, Heisenbergstmsse 1, D-70569 Stuttgart, Germany 
^SUPA, School of Physics and Astronomy, University of St. Andrews, 

North Haugh, St. Andrews, Fife, KYI6 9SS, United Kingdom 
^ Lahoratorio de Bajas Temperaturas, Departamento de Fisica de la Materia Condensada, 
Institute de Ciencia de Materiales Nicolas Cabrera, Condensed Matter Physics Center (IFIMAC), 
Facultad de Ciencias Universidad Autonoma de Madrid, 28049 Madrid, Spain 
^ Unidad Asociada de Bajas Temperaturas y Altos Campos Magneticos, 

UAM, CSIC, Cantoblanco, E-28049 Madrid, Spain 
^ Max-Planck-Institut fiir die Chemische Physik fester Stoffe, 

Nothnitzer Strasse fO, 01187 Dresden, Germany 
(Dated: December 22, 2015) 

The order parameter and pairing mechanism for superconductivity in heavy fermion compounds 
are still poorly understood. Scanning tunneling microscopy and spectroscopy at ultra-low temper¬ 
atures can yield important information about the superconducting order parameter and the gap 
structure. Here, we study the first heavy fermion superconductor, CeCu 2 Si 2 . Our data show the 
superconducting gap which is not fully formed and exhibits features that point to a multi-gap or¬ 
der parameter. Spatial mapping of the zero bias conductance in magnetic field reveals the vortex 
lattice, which allows us to unequivocally link the observed conductance gap to superconductivity 
in CeCu 2 Si 2 . The vortex lattice is found to be predominantly triangular with distortions at fields 
close to ~ 0.7 Hc2. 

PACS numbers: 74.55.+v, 74.70.Tx, 74.25.Uv 


Superconductivity in heavy fermion materials was first 
observed in CeCu 2 Si 2 [ll and was unexpected. The for¬ 
mation of heavy fermion bands is usually ascribed to 
the interaction between delocalized conduction electrons 
and localized magnetic moments, whereas localized mag¬ 
netic moments in conventional superconductors rapidly 
destroy superconductivity. Since this discovery, super¬ 
conductivity has been found in a range of other heavy 
fermion materials, often in close proximity to a quantum 
phase transition between a magnetically ordered phase 
and a phase dominated by Kondo screening. This prox¬ 
imity of superconductivity to a magnetic quantum crit¬ 
ical point found in many Ce-based compounds [H indi¬ 
cates that magnetic fluctuations and the influence of the 
quantum critical point might promote superconductiv¬ 
ity in these materials. The physics near the quantum 
critical point sensitively depends on the balance between 
screening of the local magnetic moments and interactions 
between themj^. In CeCu 2 Si 2 , slight changes in the ex¬ 
act composition (specifically the Cu-to-Si ratio) result in 
superconducting crystals {S type), magnetically ordered 
samples {A type) or samples which exhibit competing 
phases {A/S type) 0, [H. Superconductivity occurs close 
to a spin density wave-type quantum critical point and is 
expected to be unconventional in nature. Neutron scat¬ 
tering data indicate that superconductivity is mediated 
by spin fluctuations rather than by phonons @. De¬ 
spite intense research, the precise form of the supercon¬ 
ducting order parameter remains elusive. Measurements 
of the relaxation rate of the nuclear quadrupolar reso¬ 


nance 0, Q and specific heat under pressure @ have 
yielded evidence for an unconventional order parame¬ 
ter. Angle-resolved resistivity measurements of the up¬ 
per critical field Hc 2 exhibit a fourfold symmetry of Hc 2 
consistent with a dxy symmetry of the order parameter 
US- In contrast, recent thermodynamic measurements 
offer evidence for nodeless multiband superconductivity 
in CeCu 2 Si 2 , challenging the view that the pairing sym¬ 
metry is of nodal d-wave type [Uj . Also from theory, dif¬ 
ferent symmetries of the order parameters and different 
coupling mechanisms have been proposed, e.g. a d-wave 
symmetry for superconductivity mediated by magnetic 
fluctuations 13 and, very recently, 5±-wave symm etry 
emerging from magnetic fluctuations of higher order [l3|. 

Scanning tunneling microscopy and spectroscopy has 
been successfully employed to study heavy fermion 
materials [13. [l5| as well as Kondo lattice compounds [HI, 
ITj . Yet, due to the low temperatures required for these 
experiments, STM and STS have only been applied to few 
heavy fermion superconductors successfully |l8l423j. de¬ 
spite these having been thoroughly studied using macro¬ 
scopic techniques (see, e.g., ref. [2J for a review). In par¬ 
ticular, quasiparticle interference imaging provides evi¬ 
dence for a dx‘ 2 -y 2 pairing symmetry in CeCoIns [Ulii. 
Imaging the vortex lattice and transitions in the vortex 
lattice arrangement can yield additional insight into the 
structure of the superconducting gap [22| . Here we use a 
dilution refrigerator-based STM to study the supercon¬ 
ducting properties of CeCu 2 Si 2 at temperatures down to 
20 mK [ 23 . Experiments have been performed on an 
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FIG. 1. (a) Body centered tetragonal crystal structure of 

CeCu 2 Si 2 . (b) Topographic image of the surface after cleav¬ 
age showing terraces which exhibit atomic periodicities and 
step heights of ~ 2A (30x30 nm^, T = 4.5K), (c) high energy 
normalized tunneling spectrum taken at T = 200mK show¬ 
ing a W-shaped feature close to the Fermi energy (I/rms = 
3mV). (d) Normalized tunneling spectrum recorded at 20mK 
(Vrms = 40/iV), exhibiting a gap-like feature with two shoul¬ 
ders symmetric with respect to zero bias. Blue arrows indicate 
the main gap, red arrows the low energy shoulder. Red line 
shows the best fit obtained for two gaps, one of d-wave sym¬ 
metry at low energies and one of s-wave symmetry at high 
energies (for details see text and Ref. [H) . 


S type single crystal with a superconducting transition 
temperature Tc = 0.59K and an upper critical field along 
c of /io^c 2 ^ 2.3T. 

CeCu 2 Si 2 crystallizes in the ThCr 2 Si 2 M/mmm body- 
centered tetragonal structure with lattice parameters 
a = 4.1031(2)A and c = 9.9266(5)A 0, see Fig. [IJa). 
Many materials with the same crystal structure cleave 
easily perpendicular to the crystallographic c-axis (e.g. 
the 122 iron pnictides, YbRh 2 Si 2 , CeRh 2 Si 2 , YbCo 2 Si 2 
and URu 2 Si 2 ), due to differences in the bond strengths 
favoring breakage of the crystal within specific planes. 
CeCu 2 Si 2 crystals do not cleave in a way which yields 
extended, atomically flat terraces. The bonds within 
the unit cell appear to be of rather similar strength 
which prevents the development of a preferred cleav¬ 
ing plane [ 2 ^. Out of more than thirty single crystals 
which we cleaved, only two resulted in images which ex¬ 
hibit atomic periodicity and atomically flat terraces as 
reported below. 

To promote cleaving in the ab plane, a groove was cut 
around the perimeter of the sample in a plane perpendic¬ 
ular to the c axis. In order to avoid surface contamination 
the sample was cleaved in-situ in cryogenic vacuum at a 
cleaving stage fixed to the 4K plate of the cryostat. After 
cleaving, the sample was immediately transferred into the 
STM head. Differential conductance spectra were mea¬ 


sured with a lock-in amplifier, with a bias modulation 
of between 10 and 40 /iVRMS (unless stated otherwise). 
The base temperature of the instrument is below 20mK, 
the electronic temperature of the instrument has been 
determined previously to be 137mK (iBj . 

After successful cleavage of the sample, flat regions in 
an area of several hundred nanometers were found. In 
the absence of a natural cleavage plane, different cleaves 
result in different surface terminations, and also across 
the sample, the termination is not the same everywhere. 
Fig. mb) shows an example of atomically flat terraces 
with steps of a height of 2±0.2A in between. For the spe¬ 
cific area shown in fig. [11(b), we can analyze the atomic 
periodicity and tilt angle and can attribute it to a (105) 
surface. Spectra taken in a range of ±80mV show a W- 
shaped feature close to the Fermi level(see Fig. HKc)), 
reproducible across different surface terminations. It is 
suggestive to interpret the peak-like feature near zero bias 
as a signature of the heavy 4f-bands, its position is consis¬ 
tent with the energy of the Kondo resonance determined 
from angular resolved photoemission . 

At the lowest temperature achieved, 20mK, the dif¬ 
ferential tunneling conductance measured on this sur¬ 
face reveals a gap-like feature on an energy scale of 
100/reV which we attribute to the superconducting gap, 
see Fig. [11(d). The spectra do not show a strong depen¬ 
dence on the location on the surface. The tunneling spec¬ 
trum of the superconducting gap exhibits a rather com¬ 
plex structure with weak coherence peaks at ±170/iV, ad¬ 
ditional shoulders at ±75/iV and a substantial zero bias 
conductance of about 65% of the normal state differential 
conductance. We have attempted to describe the data by 
different models for the superconducting gap (for details 
of the fits see [lH). We find an excellent agreement for 
two gaps, one which describes the low energy features 
by a superconducting gap with a d-wave symmetry and 
a second one describing the high energy features by a 
gap with s-wave symmetry. The extracted fitting pa¬ 
rameters give Ai = (56 ± 3)/ieV for the d-wave gap and 
A 2 = (147 ± 4)/reV for the s-wave gap. These findings 
point towards a rather complex multi-gap superconduct¬ 
ing order parameter. The ratio of the gap magnitudes of 
the large and the small gap obtained from the best fit to 
the tunneling data A 2 /A 1 ^ 2.6 compares well with the 
ratio of 2.5 deduced from specific heat [ll|. The shape 
of the tunnelin g s pectra is similar to the ones observed 
in CeCoIn5[2i,l23|. 

Fig. Ha) shows for comparison tunneling spectra taken 
at T = 200mK (< Tc) at zero magnetic field and at 
a field of jHoH = 2T parallel to c, close to the upper 
critical field Hc 2 0. As can be seen from the spectrum 
measured at 2T, the gap has almost closed completely. 
The temperature dependence of the spectra as shown in 
Fig.[2](b) reveals that the gap vanishes as the temperature 
rises from 20mK to 600mK. 

To assess the temperature dependence of the gap size 
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FIG. 2. (a) Differential conductance spectra (normalized at 
—0.4mV recorded at 200mK in zero field (open black symbols) 
and in a field of 2T (solid red symbols). (b) Scanning tunnel¬ 
ing spectroscopy as a function of temperature (I/rms = 24/iV, 
except at 20mK Vrms = 40/iV). Red lines represent fits with 
the Dynes equation for a single s-wave gap [^; spectra have 
been shifted vertically for clarity, (c) Gap size as function 
of temperature extracted from fits in (b), the red line rep¬ 
resents the mean held behaviour expected for an s-wave or¬ 
der parameter, A(T) = Aotanh(^^^ — 1). The ht yields 
Ao = 67 ± 3/ieV and Tc = 670 ± 30mK. 


in a more quantitative way, we have used the Dynes equa¬ 
tion [ 29 I to fit the spectra in Fig. EJb) with a single 
isotropic s-wave gap. While a multi-gap order parameter 
would be more appropriate, the fit becomes unstable as 
the features of the spectra become broadened at higher 
temperatures. The broadening parameter F has been 
fixed to its value at base temperature, F = 71/ieV and 
is likely dominated by gap anisotropy and the multigap 
structure. In addition we have accounted for broadening 
due to the finite temperature and resolution. While the 
fits yield overall good agreement with the data, there are 
also clear deviations, especially for the spectra obtained 
at the lowest temperatures. We attribute both the fit 
quality and the large broadening parameter to the fact 
that the superconducting order parameter in CeCu 2 Si 2 
is more complex than just a single s-wave gap. We ex¬ 
pect the characteristic gap size extracted in this way to 
be qualitatively representative of the behaviour of the 
true gap. The temperature dependence of the gap size 
A(T) obtained in this way is shown in Fig. [2](c). It re¬ 
veals similarity to a mean-field BCS behaviour (red solid 
line) [ 3 ^. Extrapolation of the gap size to zero Kelvin 
(see fig.[2jc)) yields a value of 67±3/reV, slightly smaller 
than what is found from an isotropic s-wave fit at base 



FIG. 3. (a-e) Maps of the zero bias tunneling conductance 

(ZBTG) at 0.2 K (150x150 nm^, Erms = 40/iV) taken after 
cooling the sample in zero field as a function of magnetic field 
(0.5, 1, 1.5, 1.6 and 1.7 T), together with the autocorrelation. 
The data have been filtered for better contrast. The trian¬ 
gular vortex lattice is clearly seen, (f) Distance between the 
vortex cores as a function of magnetic field. The points are ob¬ 
tained from averaging the distances in the three high symme¬ 
try directions from the autocorrelation. Solid lines show the 
expected dependencies for a triangular (red line) and square 
lattice (blue line), (g) Anisotropy of the vortex lattice ex¬ 
tracted from the autocorrelation of differential conductance 
maps as a function of magnetic held. As a measure for the 
anisotropy, we use y (see eq. [T] inset shows the anisotropy ex- 
emplarily for the autocorrelation of the map taken at 1.6T). 
Maps taken with a tunneling set point of = 5mV (except 
(d) with U = 3mV), I = 0.6nA. 


temperature (A = 71/ieV at 20mK). The difference of 
the temperature dependence of the gap size for s-wave 
and d-wave order parameters is rather minor, so that we 
would not expect to be able to distinguish them from the 
temperature dependence of the gap sizejdlj. 

Spatial maps of the zero bias tunneling conductance 
(ZBTG) at finite magnetic fields can reveal the vortex 
lattice[32|. Observation of the vortex lattice can provide 
further information on the order parameter. For isotropic 
order parameter and electronic structure, the vortices 
usually arrange themselves in a triangular Abrikosov 
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lattice[33|. Other symmetries of the vortex lattice can 
occur if there is an anisotropic interaction between the 
vortex cores, e.g. due to the symmetry of the super¬ 
conducting order param eter or anisotropies of the elec- 

Maps of the differential con- 


M 


tronic structure 
ductance at zero bias were taken at magnetic fields of 
IIqH = 0.5,1,1.5,1.6,1.7T, see Fig. [3](a to e). As ex¬ 
pected, the number of vortices is proportional to the 
applied magnetic field. Despite the rather rough sur¬ 
face morphology, our data show an ordered vortex lattice 
which is close to triangular symmetry. 

We have analyzed the distance between the vortex 
cores ai from the autocorrelation of the ZBTC maps 
along the high symmetry directions of the vortex lattice. 
The average distance between the neighboring vortices 
as a function of magnetic field reveals a behaviour con¬ 
sistent with a triangular lattice (fig. [3](f)). As can be 
seen from figs. [3](a)-(e), our data indicate a distortion of 
the vortex lattice away from regular triangular symme¬ 
try, which becomes largest at magnetic fields of 1.6T. At 
this field, the vortex lattice also appears with a different 
orientation than at smaller or larger fields. We define 


X = max(ai, 02, a 3 )/min(ai, 02,03) ( 1 ) 

as a measure for the asymmetry of the vortex lattice, with 
X = 1 corresponding to a regular triangular vortex lattice 
(see inset in fig.[3l[g)). The magnetic field dependence of 
X is plotted in fig.ISjg), confirming a significant distortion 
of the vortex lattice specifically at a field of l.GTjs^. 
This distortion is robust across multiple data sets taken 
with different parameters, excluding drift as a possible 
reason. Also the surface roughness is unlikely the cause 
for the anisotropy, because this is similar for the images 
in fig. [H^a-e). 

Further insight into the properties of superconductiv¬ 
ity can be obtained by analyzing individual vortex cores. 
In fig. HKa), we show the decrease of the zero bias con¬ 
ductance, cr(r, 0), as a function of distance r from the 
center of the vortex core, as superconductivity recovers. 
The characteristic length scale ^ of this recovery, which 
is a measure of the superconducting coherence length, is 
determined by fittin g an exponential decay to the radial 
profile of the ZBTC [lil, 

cr(r, 0) = (Jo • -h c, (2) 

where ao is the additional conductance in the vortex core 
and c accounts for a constant background. We obtain 
^ = 7.1 ± 0.3nm averaging the decay lengths obtained 
from the data shown in fig[4|(a). This value of ^ is close 
to the one obtained previously for the coherence length 
flQnm) j39|. 

For superconductors with an isotropic gap and which 
are in the clean limit, strong vortex bound states are 
expected near the Fermi level which give rise to a 
strong peak in differential conductance spectra recorded 



FIG. 4. (a) Evolution of the azimuthal average of the nor¬ 

malized ZBTC cr(r, 0) at 0.2 K and IT, 1.5T and 1.6T as a 
function of the radial distance r from the vortex center (er¬ 
ror bars represent the standard deviation of the data, average 
over multiple vortex cores). From fits of Eq. [2] to the indi¬ 
vidual curves we obtain ^ = 7.1 ± 0.3nm on average (solid 
line), (b) Azimuthally averaged tunneling spectra as a func¬ 
tion of distance from a vortex core measured at IT. The zero 
bias tunneling conductance is substantially larger inside the 
vortex core than away from the vortex core (data has been me¬ 
dian filtered for clarity, spectra normalized at a bias voltage of 
350/iV). (c) tunneling spectra from (b) in the vortex core (red 
triangles) and away from a vortex core (black squares). Even 
in the vortex core, the zero bias tunneling conductance stays 
below the normal state value (i.e. at larger bias voltages). 


in the center of a vortex core[32j. Spectra of the differen¬ 
tial conductance acquired inside a vortex core and away 
from it show no evidence for a vortex bound state (see 
Fig.[4|(b, c)). This behaviour is consistent with a not fully 
formed gap on some bands or a considerable amount of 
scattering in the superconductor (or at the surface). The 
mean free path in CeCu 2 Si 2 is on the order of IQnmj^. 
therefore the material is not in the clean limit. 

The zero bias conductance extracted from tunneling 
spectra at zero field agrees reasonably well with specific 
heat data of Kittaka et al [Uj: at zero field and 200mK, 
our tunneling spectra yield a ZBTC of 0.74 ± 0.05 of the 
normal state conductance, whereas the normalized spe¬ 
cific heat at 200mK (Cs/Cn) is about 0.65 of the one 
in the normal state. This shows that the superconduct¬ 
ing gap we detect in tunneling spectra at the surface is 
consistent with what one would expect from bulk super¬ 
conductivity as measured by specific heat jdlj . 













5 


It is interesting to compare our observations to the case 
of CeCoIns, where experiments suggest a (i-wave sym¬ 
metry of the order parameter [23| . The spectra we 
have observed at base temperature (see fig. [Hd)) show 
a striking similarity with tunneling spectra of the su- 
perconductin g ga p in the heavy fermion superconductor 
CeCoIus [ 2 ^. 123|. Small angle neutron scattering of the 
vortex lattice in CeCoIns has previously shown substan¬ 
tial variations from triangular to square and rhombic lat¬ 
tice symmetries jiij. 

In conclusion, our tunneling spectroscopic measure¬ 
ments on the first discovered heavy fermion supercon¬ 
ductor, CeCu 2 Si 2 , show clear evidence for superconduc¬ 
tivity with at least two distinct gaps, with the best fits 
obtained for the smaller with d-wave character and the 
larger with 5 -wave character. Although the observed gap 
does not open fully, the temperature and magnetic field 
dependent measurements establish its link to supercon¬ 
ductivity. Characterization of the vortex lattice shows a 
triangular vortex lattice at low fields with significant dis¬ 
tortion at fields around 0.7idc2 indicating an anisotropic 
interaction between vortex cores. The features observed 
in the tunneling spectra at low energy support a multi¬ 
gap scenario. The shape of the gap detected in tunneling 
spectroscopy as well as the distortion of the vortex lattice 
indicate similarities with superconductivity in CeCoIns. 
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